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Abstract 



^vq In the optical, the spectrum of symbiotic binaries consists of contributions from the cool giant, symbiotic nebula and the hot 
—. star. Strong emission lines are superposed on the continuum. In this paper we introduce a simple method to extract individual 
f^j components of radiation from photometric UBV magnitudes. We applied the method to classical symbiotic stars AX Per, AG Dra, 
AG Peg and Z And, the symbiotic novae RR Tel and V1016 Cyg and the classical nova VI 974 Cyg during its nebular phase. We 
! estimated the electron temperature and emission measure of the nebula in these systems and the V magnitude of the giant in the 
^vq symbiotic objects. Our results are in a good agreement with those obtained independently by a precious modelling the UV-IR SED. 

i 1 Keywords: Stars: binaries: symbiotic; techniques: photometric 

PACS: 97.30.Eh, 97.30.Qt, 97.80.Gm 
CO 

Oh 
i 

o 

CO 

03 



> 
O 



o 
o 



X 



1. Introduction 

Symbiotic stars are long-period (P or b ~ 1 — 3 years or more) 
interacting binary systems, which comprise a late-type giant 
(usually a red giant of the spectral type M) and a hot compact 
i star, most probably a white dwarf. If the giant of the spectral 
type G or K is present, we call them 'yellow symbiotics'. Dur- 
ing quiescent phases, the white dwarf accretes a fraction of the 
wind from the giant. Typical mass loss rate from the giant in 
the symbiotic binaries is of a few times 1CT 7 M Q yr _1 . The ac- 
cretion process leads to heating up the white dwarf's surface to 
a very high temperature of 7), ~ 10 5 K, and increases its lumi- 
nosity to Lh ~ 10 2 - 10 4 L Q . Such the hot and luminous WD is 
capable of ionizing neutral wind particles of both the stars giv- 
ing rise to a strong nebular radiation. During quiescent phases, 
the symbiotic nebula represents mostly the ionized part of the 
stellar wind from the giant. During active phases, the mass loss 
rate from the hot component increases, and can temporary ex- 
ceed that from the giant. As a result, properties of the symbiotic 
nebula significantly change during the active phases. 

Accordingly, the observed flux, (corrected for the inter- 
stellar extinction) is given by the superposition of three basic 
components of radiation - from the nebula, F n 1 ehula , the cool gi- 
ant, Ff m , and the hot stellar source, F h °', i.e. 



T?nebula , j^S'int r-hot 
r A +r A +r A ■ 



(1) 



Contributions of the individual components of the radiation de- 
pend on the wavelength. The hot stellar source dominates the 
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far-f/ V spectrum, while the nebula is usually dominant in the 
neai-UV/U domain, and the cool giant (especially the red gi- 
ant) dominates the near-//? spectral region. In the case of yel- 
low symbiotics the contribution of the giant to the Johnson's 
V filter is very strong relatively to that from the nebula. How- 
ever, contributions from different sources in a symbiotic system 
depend also on the level of the activity. 

During active phases of symbiotic binaries with a high or- 
bital inclination, we observe narrow minima in their light curves. 
They can be interpreted as eclipses of the hot object by the cool 
giant. This implies that during the activity the hot stellar source 
produces a dominant amount of its radiation within the opti- 
cal region. It is believed that an optically thick disk is formed 
around t he hot star at the orbital plane during active phases (see 
Skopal, 2005). Its relatively cool rim mimics a warm pseu- 
dophotosphere, while the circumstellar material above/below 
the disk can easily be ionized by its hot central part. There- 
fore the observed depths of the minima are partially filled in by 
the nebular radiation. Thus during the eclipses, the observed 
light consists of the contributions from the giant and the neb- 
ula. During quiescent phases the optical region is dominated 
by the extensive nebula, which thus prevents from observing 
eclipses. Corresponding colour indices are very different from 
those observed during active phases. 

The aim of this paper is to propose a method of disentan- 
gling photometric U BV magnitudes of symbiotic stars into their 
components from contributing sources as given by Eq. ([TJ. 
This means to determine physical parameters of their radiation. 
In particular, the electron temperature, T e , and emission mea- 
sure, EM, of the symbiotic nebula, and the brightness of the 
cool giant in the V passband (V g,ant ). We describe our method 
in the following section. First, we formulated some simplify- 
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ing assumptions (Sect. 2.1), and then we derived a system of 
equations, which relates theoretical properties of the radiative 
components to their observational characteristics (Sect. 2.4). In 
Sect. 3 we applied this approach to selected objects and com- 
pared the results with those achieved independently by another 
methods. Discussion and summary of the results are found in 
Sects. 4 and 5, respectively. 

2. The method 



2.1. Assumptions 

In our approach we assume that the contribution from the 
hot stellar source can be neglected within the optical. We thus 
consider contributions only from the nebula and the cool giant, 
which simplifies Eq. (Q]) to F ^ - F" ehula + Ff ant . This assump- 
tion can be applied for systems during quiescent phases, while 
during active phases, it is valid only under specific conditions. 
In the following two points we summarize reasons and condi- 
tions, under which this assumption can be applied. 

(i) During quiescent phases the hot star contribution can 
be neglected in the optical, because of its very high temper- 
ature, which moves the maximum of its radiation to the ex- 
treme ultraviolet or supersoft X-ray domain, while the opti- 
cal is contributed only with a ver y faint Rayleigh- Jeans tail 
in a blackbody approximation (e.g . iKenvon & WebbinM 1 1984 ; 
Murs et et all Il99ll ISkopalL 120051) . We investigated this case 
for three classical symbiotic stars AG Dra, AG Peg and Z And 
during their quiescent phases. 

(ii) During active phases, a large optically thick disk is cre- 
ated around the accretor at the orbital plane. Its signatures 
are indicated best f or the systems with a high orbital inclina- 
tion (see Fig. 27 of ISkopali 120051) . For the eclipsing symbi- 
otics the hot stellar source is represented by a relatively warm 
(7h ~ 22 000 K) disk's rim, which light contribution can be sig- 
nificant within the optical region. Therefore, the radiation from 
the hot active object can be neglected only during the eclipses. 
As example here, we analysed the light from the 1994 eclipse 
of AX Per during its 1989-95 active phase. In the case of non- 
eclipsing systems, the optical light from the hot object can be 
neglected at any orbital phase, because we are viewing the ac- 
tive hot star more from its pole, and thus its radiation is not 
cooled by the disk material on the line of sight. We demon- 
strate this case on example of AG Dra. 

Finally, we note that our method is applicable to any object, 
whose continuous spectrum consists of two or one (nebular) 
component of radiation. We demonstrate this case on examples 
of symbiotic novae RR Tel, V1016 Cyg and a classical nova 
V1974 Cyg during its nebular phase. 



Table 1 : Colour excesses, distances and spectral types of the giants. 



Object 


Eb-v 


d 


Spectral 


ref. 




[mag] 


[kpc] 


type 




AX Per 


0.27 


1.73 


M4.5 


1,7 


AG Dra 


0.08 


1.1 


K2(K4) 


1.7 


AG Peg 


0.10 


0.80 


M3 


1,7 


Z And 


0.30 


1.5 


M4.5 


1,7 


RR Tel 


0.10 


2.5 


M6 


2,3,7 


V1016Cyg 


0.28 


2.93 


M7 


4,5,7 


VI 974 Cyg 


0.32 


1.77 
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to determine the true continuum using the multicolour photom- 
etry, we corrected the dereddened magnitudes for the influence 
of emission lines. For t he purpose of this paper we used cor- 
rections presented by Skopall l 20071) . if not specified otherwise. 
The corrections are summarized in Table [3] 

Finally, according to the Pogson's equation, we converted 
the magnitudes of the true continuum, mx, to fluxes by 



10 



-QA{m A +q x ) 



(2) 



where the constant defines the magnitude zero. For the stan- 
dard Johnson U BV photometric system and the fluxes in units 
of erg cirT 2 s~' A" 1 , qy _ 20.9 mag, c/b = 2 0.36 mag and qy = 



21.02 mag (e.g. Henden & Kaitchuck, 1982). 



2.3. Contributions from the nebula and the giant 

Now, let us have a look on basic relations for the nebular 
continuum and that from the giant. The nebular flux in Eq. (Q]) 
can be approximated by 



Ff ula = k n xs A (T e ), 



(3) 



where k n [cm -5 ] is the scaling factor, which determines amount 
of the nebula and SA(T e ) is the volume emission coefficient 
[erg cm 3 s _1 A -1 ], which depends on the electron tem perature 



2.2. The true continuum from UBV magnitudes 

Photometric measurements i.e. observed UBV magnitudes 
used in this work are summarized in Table [2] 

First, we corrected the observed magnit udes for the inter- 
stellar extinction using the extinction curve o flCardelli. Clayton & 
1989) with appropriate Eb-v colour excess (TableQ]). Second, 



of the nebula T e , and is a function of the wavelength (e.g. lBrown & M athe 
Il970l) . For the sake of simplicity, we calculated the volume 
emission coefficient for the hydrogen plasma only, including 
contributions from recombination and bremsstrahlung. In ad- 
dition, Eq. (0 requires T e , and thus the emission coefficient 
£i(7e)> to be constant throughout the nebula. The total emis- 
sion produced by the optically thin nebula is 



And l F 



2 r?nebula 



= s.i 



I n e n. 
Jv 



dV = s A EM, 



(4) 



EM [cm -3 ] is the so-called emission measure. It is deter- 
mined by the volume of the nebula and concentrations of elec- 
trons and ions (protons), n e and n + , respectively. So, with the 



2 



Table 2: Input parameters for investigated objects - Photometric measurements. 



Object 


Tnlinn rlntp 

J Lillu.ll LluLC 


u 


g 


v 




JD - 24... 


[mag] 


[mag] 


[mag] 


AX Per 


49 571 


12 689 


1 9 R7^ 


1 1 583 






19 719 
1Z. / 1Z 


1 9 SQS 






AQ SQ^ 


12 744 


12 946 


1 1 790 




49 601 


12.733 


12.936 


1 1.681 




49 608 


12 742 


12 999 


1 1 722 




49 621 


12.715 


12.870 


11.660 




S9 7^S 


i n q^7 


1 1 090 


Q 7^S 
y. 1 JJ 




53 178 


10.952 


11.044 


9.720 




52 919 


9.189 


10.022 


9.181 


AG Peg 


49 974 


9.16 


9.76 


8.56 




49 988 


9.13 


9.73 


8.55 




49 993 


9.09 


9.71 


8.55 




49 996 


9.10 


9.76 


8.58 




50012 


8.97 


9.75 


8.56 




50024 


9.08 


9.69 


8.56 


ZAnd 


45 172 


11.41 


11.72 


10.50 




45195 


11.28 


11.74 


10.64 




45 196 


11.27 


11.77 


10.64 




45 204 


11.24 


11.77 


10.67 


RR Tel 


48 066 


9.90 


11.26 


10.84 




48 067 


9.83 


11.24 


10.84 


V1016Cyg 


45 975 


10.26 


11.41 


11.15 


VI 974 Cyg 


48 883 


8.243 


9.359 


9.426 


Table 3: Input parameters - Corrections for the emission lines. 


Object 


Julian date 


AU 


AB 


AV 




JD- 24... 


[mag] 


[mag] 


[mag] 


AX Per 


49 600 


-0.77 


-0.62 


-0.19 


AGDra 


52765 


-0.033 


-0.10 


-0.01 




53172 


-0.030 


-0.08 


-0.01 




52919 


-0.081 


-0.12 


-0.02 


AG Peg 


49305 


-0.19 


-0.41 


-0.12 


ZAnd 


45 293 


-0.30 


-0.48 


-0.12 


RR Tel 


51 836 


-1.71 


-1.68 


-0.87 


V1016Cyg 


46 045 


-1.23 


-1.67 


-1.19 


VI 974 Cyg 


48 883 


-1.57 


-1.73 


-0.66 



aid of Eq. (13), the emission measure can be expressed as 



To simplify relations introduced in the following section, we de- 
note the colour indices of the giant's spectrum, U-B and B - V, 
by UB and BV, respectively. 

2.4. Parameters of the nebula and giant 

Here we propose a method to estimate parameters of the 
nebular radiation (T e , EM) and that of the giant (e.g. V gUmt ), 
which dominates the optical part of the spectrum of most sym- 
biotic stars. Based on the assumptions formulated in Sect. 2.1. 
and relations for the nebular emission (Sect. 2.3.), we can write 
following set of equations for the considered components of ra- 
diation. First, the total flux in the continuum is given by the 
superposition of fluxes from the nebula and the cool giant, i.e. 



r-nebula r<giant 



10 



-0.4(U c ""+q v ) 



(7) 



r^nebula , pgiant 



ehula 



^Q-Q.4(B c ' m, +q B ) 
10 -0.4(V"'"'+ 9v ) 



(8) 

v + Ff" = io- - 4 ^"'"^), ( 9) 

where U con ', B cont and V cont are magnitudes of the true contin- 
uum (Sect. 2.2.). Second, according to Eq. (0, the nebular flux 
in the spectrum is a function of the scaling factor k n and the 
electron temperature T e , i.e. 

^nebula 



U 

nebula 
B 

nebula 
V 



■ k n s B (T e ), 
k n £\r(T e ). 



(10) 

(11) 
(12) 



Third, the ratio of the fluxes from the giant in different filters 
can be expressed with the aid of the Pogson's equation as 



7 gitmt 



r^giant 



10 



-0.4(f/ s " , "'-B' !io "'+ 9 u-?B) 



10 



-0.4(UB+q v -q B ) 



pgiant 

F B _ l Q-QA{B«"-"'-V^"'+q & -q w ) = , n -0.4(BV+ ?B -?v) 



7 gianl 



IQ- 



(13) 



(14) 



Thus we have a system of eight equations (Eqs. (0 to ( fl4b ) for 
eight variables, 



^nebula 



^nebula r'nebula 



pgiant 



pgianl 



pgianr 



k n and T e 



^nebula 



EM = 4nd 2 k n = And 



2 A 



(5) 



which can be solved for the known spectral type of the cool 
giant (i.e. the U-B and B - V indices) and the measured 
UBV magnitudes. The aim is to determine the three fundamen- 
tal parameters, the electron temperature T e , the scaling factor 
k n and the giant's magnitude, V s "" 1 ' . The first two parameters 
determine the nebular emission, while the third one settles the 
radiation from the giant characterized with indices U-B and 
from both sides of the Balmer jump rsee lSkopalLl2009h. For the B - V, and thus defines the second term on the right side of Eq. 
respon se function of the U filter, as published bv lMatthews & Sandafl t- From the abovementioned system of equations (Eqs. © 

to (fl4l ) we can derive a relation for determining the electron 
temperature in a form (see Appendix A) 



where d is the distance to the object. 

Due to the asymmetry of the U filter with respect to the 
wavelength of the Balmer discontinuity, /i-Baimer, we calculated 
the emission coefficient s\j as the weighted avera ge of i ts values 



(1963), we get 



£U = 0.6 £u- + 0.4 £(j+ 



(6) 



where E\j- and £u+ are emission coefficients at the short and 
long wavelength side of the /iBaimer, respectively. 

The stellar radiation from cool giants of different spectral 
types can be characterized by different colour i ndices. For the 
purpose of this work we used colour indices from ljohnson ( 1966b . 



£ B(T e ) r q-0.4(V"""+£/B-?b 



10 



-0.4(U c ""'-BV-q B )] 



[ 10 -0.4(E/' 



10 -0.4(V"»"- 9u )l 



^q-0.4(B"""-BV- 9u ) 

9v) _ iQ-0A(B c ""'+UB-qy)j _ q 



(15) 



3 



Solving this equation for T e , allows us to determine easily other 
parameters. For example, the scaling factor k a as 



pnebula ^Q-0A(BV+V"""+q B ) _ iq-0 A(B"""+q B ) 



s Y (T e ) io-°- 4 < Bl/+ ?B-?v) fiv (r e ) - £B (r e ) 

and the giant's magnitude in the V passband V imnt as 



(16) 



ygiant = _ 2 5 log ^ 10 -0.4(V»"' +9V ) _ pnebula j _ ^ , , ? , 



where Fy bula is given by 



^nebula 



lQ-0A(BV+V< m "+q B ) _ iq-0.4(B™"+ 9b ) 
lQ-OA(BV+q B -qv) _ £ B( r c) 



(18) 



For more details see Appendix A. The ratios of the emission 
coefficients eu/sv and sb/ev are plotted in Fig.Q] 
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Figure 1: Ratios of the emission coefficients we used to solve Eq. 4151 . 



3. Application to selected symbiotic stars 

Our resulting parameters, the electron temperature T e (from 
Eq. (fl3Tl), scaling factor k a (Eqs. Q and (fToll). emission mea- 
sure EM (see Eq. (0) and the V magnitude of the giant V gumt 
(Eq. ( TT7i i) are summarized in Table [4] Values of EM were 
scaled for distances summarized in TableQ] Disentangled UBV 
magnitudes for the investigated objects are plotted in Fig. 2. 

Assuming that the uncertainties of the observed U, B and V 
magnitudes for all investigated syst ems are ±0-05 mag, +0.03 mag 
and +0.03 mag, respectively (e.g. ISkopal et all 120071) . we de- 
rived intervals of possible values of the searched parameters, 
which can be found in the fourth, sixth, eighth and tenth col- 
umn in Table |4] Uncertainties in the corrections for emission 
lines are not included. If more photometric measurements were 
available, we used mean values of parameters, corrensponding 
to the individual UBV magnitudes. For these cases, we adopted 
uncertainties as the interception of parameter ranges given by 
individual measurements (except for the brightness of the cool 



giant in Z And and AX Per where we took the mean values of 
upper and lower limit of intervals because those interceptions 
were empty). Investigated objects are discussed in detail in the 
following subsections. 

3.1. Classical symbiotic stars during quiescent phase 
3.1.1. AGDra 

AG Dra is a yellow symbiotic star, beca use its cool compo- 
nent w as classified as a normal K2 giant (e.g. Mii rset & Sch mid. 
1999h . To analyse a quiescent phase of AG Dra by o ur method, 
we us ed U, B, V measu rements made on 200 3/05/05 dSkopal et al 



2004) and 2004/06/21 (ISkopal et a l„ 2007). We selected these 



dates, because of having simultaneous spectroscopic observa- 
tions, which allowed us to determine corrections for emission 
lines on the UBV magnitudes (Fig. 6 of lGonzalez-Riestra et al. , 
20081 and Viotti, private communication). The colo ur indices 
for a K2 giant, B — V = 1.16 and U - B = 1.16 djohnsonl 



1966b are, however, in conflict with the observed indices of true 
continuum, (B - V) b s -1.3 and (U - fi)obs - -0.2, because 
observations contain contributions from both the giant and the 
nebula. Therefore, for the purpose of our analysis we adopted 
the giant's spectral type of K4, which is characterized with in- 
dices B - V = 1.41 and U - B - 1.66 Jjohnsonl[l966l) . 

Applying our method (Eqs. (T5[ to ( fTTT i) to these data we 
obtained the electron temperature of the nebula, T e a 20 000 K 
and its emission measure, EM ~ 1.7 X 10 59 cm' 3 . 

3.1.2. AG Peg 

AG Peg is known as the slowest symbiotic nova. Currently, 
it displays all signatures of a classical symbiotic star in a quies - 
cent phase dKenvon et al. , 1993 ; Miirset & Nussbaumei , 1994 ). 
We investigated AG Peg around its optical maximum in 1993 
November, at the orbital phase cp = 0.63. Selection of these 
photometric data was important to compare our results with 
another analys is and to use appropriate corrections for emis- 
sion lines (see ISkonall l2007h The U, B, V magnitudes were 



taken from Tomov & Tomova ( 1998). We used values of 6 mea- 



surements obtained between JD 2 449 974.3 and JD 2 450 024.3 
and then averaged our results obtained from individual nights. 
For the spectral type of the giant in AG Peg we adopted M3 
dMiirset & Schmidlll999l) . 

Analyzing these data we revealed T e ~ 17910K and EM ~ 
4.85 x 10 59 cm ' 3 , which are quite similar to those obtained by 
ISkopallri2005l ). 



3.1.3. ZAnd 

Z And is a prototype of the class of symbiotic stars. We 
analysed U, B, V measurements taken around a light maximum 
of a quiescent phase (i.e. ip ~ 0.5), because the corrections 
for emission lines were made at the same orbital position of the 
binary. For the purpose of our an alysis, we used U, B, V magni- 
tudes taken bv lBelvakina ( 1992 ). We selected 4 measurements 
from JD ~ 2 445 1 7 1 .5 1 to JD ~ 2 445 204.42. The spect ral type 
of the giant in Z And is M4.5 dMiirset & Schmidi Il999h . How- 
ever, the giant's colour indices are available only for spectral 
types of M4 and/or M5. Therefore, we used indices for a M4 



4 



and M5 giant separately, and then averaged our results obtained 
from individual nights. 

Our resulting temper ature, T e ~ 29 300 K, is somewhat higher 



than that determined by ISkopal (2005). Using indices of a M5 
giant led to a lower temperature, T e ~ 28 390 K, while indices 
for a M4 giant yielded T e ~ 30 220 K. Our emission measure, 
EM ~ 7.8 x 10 59 cm' 3 , i s quite similar to that determined inde- 
pendently dSkopall[2005h . 



3.2. Classical symbiotic stars during active phase 
3.2.1. AX Per 

AX Per is known as e clipsing symbi otic binary with an 
orbital period of 680 days dSkopali 1 199 lb . The cool compo- 
nent of the binary i s a n ormal giant of the spectral type M4.5 
(Miir set & Schmidl 1 19991) . During active phases we observe 
narrow minima in the light curve at the position of the inferior 
conjunction of the giant. It is believed that they are caused by 
eclipses of the hot object by the cool giant. 

In this paper we studied the 1994 eclipse, observed dur- 
ing the 1989-95 active phase of AX Per. We used 6 UBV 
measurements fr om the totality, betw een JD 2449 570.54 and 
JD 2449620.56 dSkopal et al.Lll995h . We assumed that contri- 



butions from both the nebula and the giant had not been varying 
significantly during the relatively short time of the total eclipse 
(~ 2 months). Therefore, also in this case, we derived parame- 
ters for each night set of U BV magnitudes, and then used then- 
means, to get the most probable values of the fitting parameters. 
As in the case of Z And we applied our method to both spectral 
types of a giant, M4 and M5, and used their means to obtain 
resulting parameters. 

In this way we received the electron temperature of the neb- 
ula J e ~ 30 770 K, the scaling factor, k n ~ 0.41 x 10 15 cnT 5 
which corresponds to the emission measure, EM ~ 1.47 x 
10 59 cm" 3 . Using the giant of the spectral type M5 led t o lower 



values of both nebular parameters than for a M4 giant. ISkopal 
J2005I) modelled the UV-IR SED from the 1990 eclipse and 
found very similar value of T e , but a factor of ~2 higher EM 
than we obtained for the 1994 eclipse. The lower amount of the 
nebular emission indicated during the 1994 eclipse is probably 
connected with the ending of the 1989-95 active phase. Note 
that the high temperature nebula appe ars to be stron g just dur- 
ing active phases of symbiotic stars (Skopal, 2005). Finally, 
the mean value of the giant's V magnitude is 11.15 mag. How- 
ever, its value runs from 1 1 .02 to 1 1 .30 mag for these 6 individ- 
ual UBV measurements, including their uncertainties. Different 
values of the V s '""' magnitudes are beyond the measured uncer- 
tainties. Therefore, we conclude that they can be caused by a 
variable brightness of the giant during the eclipse, perhaps due 
to its pulsations. 

3.2.2. AGDra 

Here we demonstrate our method on the UBV magnitudes 
measured duri ng the optical burs t of AG Dra in 2003 October 
(see Fig. 8 in lSkopal et ail 120071) . We derived T e « 45 000 K 
and£M ~ 1.9 x 10 60 cm , which are significantly higher than 
those we indicated during quiescent phase. To check sensitiv- 
ity of these parameters to the contribution from the giant, we 



used different colour indices of the giant. By this way we ob- 
tained a lower T e ~ 30 000 K and 38 000 K for a K2 and K3 
giant, respectively. However, the emission measure was rather 
insensitive to such a small change in the spectral type of the 
giant, because of its large quantity. We note that during active 
phases the emission measure of the symbiotic nebula in AG Dra 
increased by a factor of ~10. 

In both cases (quiescence and activity) our parameters were 
in a good agreement with those obtained independently by a 
precise m odelling the (0 .12 - 5) fim SED in the continuum (see 
Fig. 13 of lSkonaiEool) . 



3.3. Symbiotic novae 

3.3.1. RR Tel 

RR Tel is a slow symbiotic nova that underwen t an outburst 

in 1944 and whose optical l i ght cu rve is still fading dMiirset & Nussbaumi 
ll994tlKotnik-Karuza et al.l 12006). The binary contains a Mira- 

type variable as the cool component, which pro duces a massive 

wind and give rise to a strong dust emission (e.g ljurkic & Kotnik-Karuzal 
2007| lAngeloni et all l2007[). The system contains a ve ry hot 
and luminous white dwarf dMiirset & N ussbaumeri ll994h . which 
ionizes a fraction of the giant's wind re sulting in a strong nebu - 
lar emission in the spectrum of RR Tel dBrvan & Kwokul99ll) . 

For the spectral type o f the giant in RR Tel we adopted M6 
dMiirset & Schmidl 1 19991) . We applied our method of disen- 
tangling the UBV magnitudes to the mea surements made on 
1990/06/23 and 1990/06/24 as published bv lMunari et ail d 19921) . 
We found the electron temperature and the emission measure 
that are typical for the quiescent phases of other symbiotic stars 
(Table©. 

3.3.2. V1016Cyg 



V 1016 Cyg underwent its nova-like outburst in 1964 (Fitz gerald et al. 

19661) . The binary contains a Mira vari able as the cool compo- 
nent with a strong IR dust emission (e.g.lTaranova and Shenavrini 



2000tlParimucha et allEoollArchipova et al.Ll2008l). The spec 
tral ty pe of the gi ant in V1016 Cy g is M7 dMiirset & Schmidl 
1999). However, Johnson! (|1966) doesn't provide the indices 
for this spectral type. Therefore, firstly we used the spectral 
type M6, and secondly, we neglected the giant contribution to 
UBV magnitudes, at all. We found that there is no significant 
difference between these two possibilities. As a giant of the 
spectral type M7 is even cooler than that of the M6 type (i.e. 
its contribution to the optical is lower), we decided to neglect 
the giant contribution and calculated only the nebular contin- 
uum. Also, it is important to note that the nebular contribu- 
tion dominates the optical (see Fig. 5 of |Skopal, 2007), which 
makes it difficult to extract the relatively very faint light from 
the giant by our metho d. In our analysis we u sed UBV mea- 
surements published bv lParimucha et al. (2000) from 1984 Oc- 
tober, which are close to dates of the spectroscopic observa- 
tions used for emission lines corrections. The large emission 
measure, EM ~ 10 60 cirr 3 , is typical for the active symbiotic 
systems. 
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Figure 2: The UBV SED for investigated symbiotic binaries and novae. The UBV fluxes in the continuum, Fjj°BV (S ect - 2.2), and their disentangled components 
from nebula, F™^j a , and the giant, F^™, are denoted by symbols in the legend. The nebular continuum, F™ bu,a , is drawn with a solid line. 
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Table 4: Physical parameters of the nebula and the giant for selected symbiotic stars and novae obtained by our method of disentangling their UBV magnitudes. 



Object 


Julian date 


T e 


T ejnin ~ 


Te.max 


k n 




"■n,ma.\ 


EM 


EM mi „ 


- KM max 


ygUmt 


ygumt _ 

min 


v gtanl 
' max 


ph.* 




JD- 24... 


[K] 


[K] 


[10 15 cm" 5 ] 


[10 15 


cm~ 5 ] 


[10 59 cm" 3 ] 


[10 59 


cm -3 ] 


[mag] 


[mag] 




AX Per 


49 600 


30770 


27 600 - 


37 600 


0.41 


0.39 


-0.47 


1.47 


1.40 


- 1.68 


11.15 


11.08 - 


11.20 


A,E 


AGDra 


52765 


20710 


12 330- 


35410 


1.24 


0.70 


- 1.90 


1.80 


1.02 


-2.75 


9.59 


9.53 - 


9.66 


Q 




53 172 


19230 


1 1 200 - 


32610 


1.16 


0.62 


- 1.79 


1.68 


0.90 


-2.59 


9.57 


9.50- 


9.64 


Q 




52919 


44 940 


34 210- 


64160 


13.2 


11.2 


- 16.1 


19.1 


16.2 


-23.3 


9.49 


9.39- 


9.61 


A 


AG Peg 


49 305 


17910 


16 500- 


19130 


6.33 


6.00 


-6.93 


4.85 


4.60 


- 5.30 


8.56 


8.51 - 


8.60 


Q 


ZAnd 


45 293 


29 300 


24 370 - 


38 040 


2.89 


2.59 


-3.14 


7.78 


6.99 


- 8.45 


10.08 


10.01 - 


10.16 


Q 


RR Tel 


51836 


18910 


16 370- 


21860 


0.81 


0.73 


-0.88 


6.02 


5.44 


-6.60 


11.83 


11.74- 


11.90 


Q 


V1016Cyg 


46045 


24480 


22 100 - 


26 530 


2.59 


2.52 


-2.71 


26.6 


25.8 


-27.9 








Q 


VI 974 Cyg 


48 883 


35 760 


31480- 


41 420 


20.4 


19.3 


-21.8 


76.5 


72.3 


- 81.8 








Q 



* phase: active phase (A), eclipse (E), quiescent phase (Q) 



3.4. Classical nova V1974 Cyg 

VI 974 Cyg (Nova Cygni 1992) be longs to the gr oup of the 
classical novae. It was discovered by IColhnsI (11992) on 1992, 
February 19. To investigate the physical conditions in this nova 
by our method, we used the data from its nebular phase taken 
around day 210 after the optical maximum. A ccording to the 
X-ray light curve of V1974 Cyg (see Fig. 1 of Rra utter et al. 



ray 

m. 



1996), the burning white dwarf in the nova was already very 



hot at that time. Therefore we could neglect its contribution to 
the optical as well as that from the red dwarf companion. The 
optical spectrum was thus strongly dominated by the nebular 
radiation. As a result, we derived parameters characterizing the 
nebular radiation component only. For this purpose, we used 
two flux-points corresponding, for example, to U and B mag- 
nitudes. Particularly, we used photometric meas urements taken 



at JD 2448 883.4 (1992/09/17) as published bv lChochol eta! 
dl993h . 



Disentangling these U and B measurements, we found a rel- 
atively high electron temperature of T e ~ 36 000 K, and also a 
very high emission measure, EM ~ 8 x 10 60 cm -3 (Table 4). 



4. Discussion 

Physical parameters for the symbiotic nebulae, which we 
derived by disentangling the UBV measurements (Sect. 3, Ta- 
ble 4), are in a good agreement w ith those deter mined by a pre- 
cise modelling the UV-IR SED dSkopali [2005h . We found that 
the emission measure can be determined with a higher accuracy 
than the electron temperature of the nebula. The former is given 
by scaling the emission coefficient, and thus can be estimated 
within the uncertainty of the photometric measurements, while 
the latter is given by the profile of the volume emission coef- 
ficient. In addition, for lower values of electron temperatures, 
to say r e < 25 000 K, the profile is relatively steeper than for 
higher values of T e . As a result a lower T e can be determined 
with a better accuracy than a higher one, and vice versa. Tak- 
ing into account difference in T e between quiescent and active 
phases, we can, in general, conclude that during the quiescence, 
T e can be determined more precisely than in the activity, in spite 
of equal uncertainties of UBV measurements. 



Parameters of the nebula in our model also depend on colour 
indices of the cool giant, which values, however, are often dif- 
ferent in different artic les. For the purpose of this work we 
adopted indices from Ijohnsonl dl966h . Further, it is important 
to know the correct spectral type of the giant. Also here differ- 
ent authors recommend slightly different spectral types for the 
giants in the investigated systems. Furthermore, our analysis 
revealed that the brightness of the cool giants can vary within a 
few times 0.1 mag (e.g. AX Per during the 1994 eclipse). This 
type of variability can be ascribed to pulsations of late-type gi- 
ants, because they represent their typical behaviour. All these 
uncertainties cannot be included in our analysis. Therefore, our 
resulting magnitudes of the giant can be considered only as es- 
timates of their mean values. 

It is of importance to note that we modelled just the opti- 
cal continuum. So it is inevitable to know corrections for the 
influence of emission lines on the UBV magnitudes. Uncer- 
tainty (ignorance) of these corrections can represent a signifi- 
cant source of errors of the fitting parameters. 

Knowing actual colour indices of the giant and corrections 
for emission lines, our system of equations allowed us to study 
the effect of uncertainties of the UBV magnitudes to the fit- 
ting parameters. We could easily found extreme values of the 
parameters, as well as their values based on all possible combi- 
nations of UBV magnitudes. 

During quiescent phases we found the electron temperature 
T e ~ 20 000 K, while during active phases our method indicated 
higher values of T e ~ 30 000 - 40 000 K. During quiescence the 
emission measure was in the order of ~ 10 59 cm 4 , while during 
activity we identified its increase to ~ 10 60 cm 4 . For symbiotic 
novae we obtained even higher values of ~ 10 60 - 10 61 cm . 

Finally, we note that our method of disentangling the UBV 
magnitudes can be applied to any spectrum composed from a 
nebular and stellar component of radiation. 

5. Conclusion 

In this contribution we employed a method of disentangling 
the composite spectrum of the optical continuum on the basis 
of simple multicolour photometric measurements. Our method 
(Sect. 2) allowed us to determine the physical parameters of the 
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main contributing components of the radiation into the optical 
region - the nebula and the giant. Main results may be summa- 
rized as follows. 

(i) On the basis of the UBV photometric measurements, we 
determined the true optical continuum of the selected sym- 
biotic (-like) binaries. 

(ii) We compared the observed continuum by a model, which 
includes contributions from the nebula and the giant (Sect. 3). 
In this way we determined the electron temperature and 
emission measure of the nebula, and the V magnitude of 
the giant. 

(iii) Our model parameters are well comparable with those de- 
termined independently by another method. In particular, 
by a precise m odelling the UV-IR SED as introduced by 



Skonall(l2005h 



(iv) Our approach thus provides a good estimate of the phys- 
ical parameters of contributing sources of radiation into 
the optical on the basis of a simple U BV photometry. 



We have 7 equations now (IA^IAT21IaT41Ia31Ia1^IAT71IaT91 >. 
Further, we express the flux from the giant in the B passband 
F g ' m ' from Eq. (K3\ as 



pgiant _ jq-0.4(B'°'"+ 9 b) _ p 



nebula 



and substitute it to remaining equations. Equations ( 1A.U . ( 1A.4I >. 
( IA.5b . (IA.6I ) don't change. The Eq. iA.71 turns into 



pgiant _ jq-0 A(UB+B ca '"+q v ) _ j Q-0.4(£/fl+(fu-gB) p 



nebula 



The Eq. JA.91 > turns into 



10 



-0.4(B"»"+ 9B ) 



t nebula 



_ |Q-0.4(BV+V a ""+4 B ) _ iQ-0A(BV+q B -qv)^. 



nebula 



(A. 10) 



(A.ll) 



We have 6 equations now (IA^IAT1Ia31|A.6IIA.10||A~TTT i. 

Further, we express the flux from the giant in the U passband 



from Eq. ( IA.U as 
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pgiant _ jQ-0.4(t/"""+9u) _ p^ 



nebula 



and substitute it to remaining equations. Equations ( 1A.4I >, dA.5b , 
( IA.6b . (1A.1U don't change. The Eq. ( 1A.10I ) turns into 



Appendix A. Relation for determining the electron temper- 
ature 

In this appendix we derive relation for determining the elec- 
tron temperature (Eq. [T5l ) from the following system of equa- 
tions 



pgiant 
U 

pgiant 
pgiant 



pgmru p, 



pgturu p, 



pg,U,U p, 

v v 

-rpnebuia 



nebula j q-0.4(J/ co "' +qv) 

nebula _ jq-0.4(S'""'+(J B ) 
nebula |q-0.4(V cwo '+^v) 

K £u(T e ), 



^nebula 
B 

^nebula 



pneouta = fe 



7 giant 



7 giant 



■pgiant 



pgiant 



= k n By(T e ), 
1 Q-0.4(Bs / '""-VS"""+ ?B -?v) _ iQ-OAW+qv-qv) 



(A.l) 

(A.2) 
(A.3) 
(A.4) 
(A.5) 
(A.6) 

(A.7) 
(A.8) 



At the beginning we have 8 equations ( 1A.1I - IAT8I ). 

We express the flux from the giant in the V passband F^""' from 

Eq. (1A.3E as 



pgiant _ jq-0.4(V"""+ 9v ) _ p^ 



nebula 



and substitute it to the remaining equations. Equations ( IA.U . 
SKU . SKM . dA3l ( TA31 i. dA^b do not change. The Eq. dA78l > 
turns into 



pgiant _ jq-().4(BV+V"""+ij b ) 



jq-0.4(W+4 B -?v) pnebula ^ 



^Q-0.4(J/™"'+(?u) _ pnebula 



10 



u 

-0A(UB+B am, +q v ) 



10 



-0A(UB+q v -q B ) jpnebula 



(A.12) 



We have 5 equations now (IAT41Ia31Ia31IA~TT1IA~T21 ). 
Further, we express the scaling factor of the nebula k n from Eq. 
( lA.4t as 



^nebula 



and substitute it to remaining equations. Equations ( IA.1 IK ( IA.12I ) 
don't change. The Eq. (lA.5t turns into 



pnebula 



The Eq. (IA.61 > turns into 



ebula 



g B(r e ) neBu i a 

s v (T e ) u • 



Ey(T e ) ne b u l a 



(A.13) 



(A. 14) 



We have 4 equations now (IA3T1lA^lA~T3llA~T4l i. 

Further, we express the flux from the nebula in the U passband 

pnebula fmm Eq (Qyj) afJ 



^nebula 



e B (r e ) J 



-nebula 



and substitute it to remaining equations. Equation dA. 1 11 doesn't 
change. The Eq. dA. 14b turns into 



^nebula 



e B (r e ) J 



-nebula 



(A.15) 



The Eq. ( IA. 12b turns into 



10 -0.4((/'""'+ ? u) _ 



Bu(T e ) 



?nebula 



£B(T e ) (A. 16) 

_ 2Q-0.4(f/B+B™"+?u) _ iQ-O.MUB+qu-qB) pnebula 

We have 3 equations now dA.llllA~T31IA.16b . 

Further, we express the flux from the nebula in the B passband 

F Ma from Eq ( T-j-g , as 



^nebula 
B 



gB(~e) 



^nebula 
V 



and substitute it to remaining equations. The Eq. ( IA. lib turns 
into 



yj-0A(B«""+q B ) _ 



-nebula 



ev(T e ) v 

2q-0.4(BV+V i °"'+(Jb) _ JQ-0.4(BV+<y B -<?v) pnebula 



(A. 17) 



The Eq. ( IA. 16b turns into 
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r 0.MU"""+qv) 



Su(~e) pnebula 

s v (T e ) v 



2Q-0.4(f/B+B ra '"+?u) _ iQ-0.4(t/fi+t?u-t? B ) g B(^e) pnebula 

sv(T e ) v 



(A.18) 



We have 2 equations now (IA. 171 [A~T~~) . 

Further, we eliminate the flux from the nebula in the V passband 

pnebula from Eqs ft£JJ\ ^ ""-gj Wg ggt the f H owing 



V 

equation: 



■^nebula _ 



jQ-0.4(fiV+V™"'+ ?B ) _ IQ-0 .4(B™"'+ 9B ) 

2fl-0.4(By+ 9B -9v) _ £ B( r e) 
«v<Te) 

10 -0.4([/B+B™"+ 9u ) _ ^q-O^CC/^'+^u) 
in-0.4(£/S+ 9u - 9B )£Bac) _ eu(r c ) 



(A. 19) 



When we multiply Eq. (IA. 19b by its denominators and do some 
other simple mathematical operations we get 

g B<Te) r^Q-0.4(£/B+BV+y™"'+ 9 u) _ 1 Q-0.4((/"""+ ?u )~| + 
Sy(T e ) L J 

+ gu(Te) r 10 -0.4(B™"+ 9B ) _ 10 -0.4(BV+y-'»'+ 9B )] + 

e v (r e ) L " J 

+ 1Q- - 4 ( bv+!/ '°"'+9li+9b-9v) _ 2Q-0.4(£/B+BV+B I ""'+c/u+?B-9v) _ Q 

For the sake of simplicity we divide previous equation by 
lQ-0A(BV+qu+q B ) _ Finally we get the equation for determining 
the electron temperature in a form 



S s(Te) r 1 Q-0.4(V'°"'+£/B-, ?B ) _ ,Q-0.M.U M -BV-qB)] , 

£v(~e) L J 



g u( r e) r 10 -0.4(B™"'-BV- ?u ) _ jq-O.^V"'"-^)] + 
+ ^Q-0.4(C/"'"'- 9v ) _ 10 -0.4(B'""'+£/S-9v)j _ Q; 



(A.20) 



which we used in this work (Eq. (15) in the main text). After 
determining the temperature r e from Eq. dA.20b . we can easily 
determine other parameters, F™ bula , F™ bula , F™ bula , k n , 



pgiant 



^gianl 



and F g ™'. 



For example, we can use equations in the 
frames. If we are interested only in three fundamental parame- 
ters (T e , k n and V 8 '""') it is better to determine the scaling factor 
of the nebula, k„, from 



pnebula 



sv(T e ) 



(A.21) 



as we did in the main text. 
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